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Abstract 

We study the properties of a parity- and time-reversal- (PT) symmetric tight-binding chain of 
size N with position-dependent hopping amplitude. In contrast to the fragile PT-symmetric phase 

g 

of a chain with constant hopping and imaginary impurity potentials, we show that, under very 

Hi 

general conditions, our model is always in the PT-symmetric phase. We numerically obtain the 
energy spectrum and the density of states of such a chain, and show that they are widely tunable. 

> 

By studying the size-dependence of inverse participation ratios, we show that although the chain 

t> 

is not translationally invariant, most of its eigenstates are extended. Our results indicate that 
tight-binding models with non-Hermitian, PT-symmetric hopping have a robust PT-symmetric 



> 

X 



phase and rich dynamics which may be explored in coupled waveguides. 
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Introduction: Since the seminal paper by Bender et al. [T] a decade ago, it has become 
clear that non-Hermitian Hamiltonians with parity and time-reversal (PT)-symmetry can 
have purely real spectra [21 E] and, with an appropriately redefined inner-product, they lead 
to orthogonal eigenvectors [2], unitary scattering jl] and, therefore, a consistent quantum 
theory. The theoretical work on continuum, PT-symmetric, non-Hermitian Hamiltonians [5] 
since then has been accompanied, most recently, by experiments in optics where spontaneous 
PT-symmetry breaking in a classical system has been observed in waveguides with PT- 
symmetric complex refractive index [HI Ej and by theoretical studies of distributed-feedback 
optical structures that can be mapped onto a relativistic, PT-symmetric Hamiltonian [Sj. 

Idealized lattice models have been popular in physics due to their analytical and nu- 
merical tractability, the absence of divergences [9], the availability of exact solutions [1QJ, 
and the ability to capture counter-intuitive physical phenomena [llj. As with the standard 
quantum theory, these models have been based on Hermitian Hamiltonians. In recent years, 
tight-binding models with a Hermitian hopping and PT-symmetric, complex, on-site poten- 
tials [TSHT^] . non-Hermitian transitions [J5], and PT-symmetric spin-chains [TB] have been 
extensively explored. For a tight-binding chain with PT-symmetric impurity potentials, a 
salient result is that its PT-symmetric phase - the range of model parameters that lead to a 
real spectrum - is extremely fragile [T2l ITi] . This fragile nature of the PT-symmetric phase 
precludes effects such as the Anderson localization [17], impurity-bound states [18] . and the 
Luttinger-liquid behavior [TH] in such a chain with a non-Hermitian Hamiltonian. 

In this paper, we explore the properties of a tight-binding chain of size A with PT- 
symmetric, non-Hermitian, position-dependent hopping amplitudes. Our main results are 
as follows: (i) We show that the system is always in the PT -symmetric phase under very 
general criteria that we derive, (ii) The energy spectrum and the resulting density of states 
in such a chain are widely tunable and symmetric around zero, (iii) Although the chain is not 
translationally invariant, (a majority of) its eigenfunctions are delocalized. Our results show 
that a robust PT-symmetric chain has non-Hermitian hopping amplitudes and Hermitian 
potentials, and that its Hamiltonian is similar to that of a chain with position-dependent, 
parity-symmetric hopping [20J. 

Tight-binding Model: We start with a Hamiltonian for an A-site tight-binding chain, 
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where c^(c„) is the creation (annihilation) operator at site n, U are the position-dependent 
hopping amplitudes, and the asterisk denotes complex conjugation. The parity operator on 
the chain is given by (m\V\n) = 5 m ^ + i^ n = 5 m>n where \m) represents a single-particle state 
localized at site m, and m = N + 1 — mis the reflection-counterpart of site m; it follows that 
Hpt, although not Hermitian, is PT-symmetric. We consider only the single-particle sector 
and, since periodic boundary conditions are incompatible with the PT-symmetry, use open 
boundary conditions. Numerical results indicate that the spectrum of Hp? is purely real 
when the hopping elements have the same sign; in the following paragraph, we analytically 
derive the criteria that guarantee this robustness. 

Let us consider a similarity transformation |5] of the non- Hermitian Hamiltonian, Hp? — > 
H = M~ l H PT M where M = diag( mi, . . . , mjv) is a diagonal matrix. It is straightforward 
to show that the transformed matrix H is Hermitian, H = H\ if and only if 

m fc+l m fc+l = ( tN-k \ > Q / 2 x 

m* k m k \ h J 

We note that this constraint only applies to a diagonal M. Thus, the PT-symmetric Hamil- 
tonian Hpt is similar to a Hermitian Hamiltonian H if and only if the phases of the hop- 
ping amplitudes (tk,tN-k) are the same for all k = {1, . . . , N — 1}. When the hopping 
elements are real, it implies that tk and t^-k must have the same sign; when they are 
complex, t m = \t m \ exp(i8 m ), it implies that 9%. = 0N-k- The eigenvalue spectrum of the 
non-Hermitian Hamiltonian Hpx is purely real, as long as these general requirements are 
satisfied. Since H = M~ x HppM = W, it follows that the eigenvalues E n of H and Hpp 
are the same, and that the orthogonal eigenvectors of H, H\v n ) = E n \v n ), and the (non- 
orthogonal) eigenvectors of the PT-symmetric Hamiltonian, HpT\u n ) = E n \u n ) are related 
by \u n ) = M\v n ). This relation provides the requisite inner-product under which the eigen- 
vectors \u n ) of Hpt are orthonormalized. We note this transformation corresponds to the 
positive-definite, self-adjoint, invertible metric i]^ 1 = MM^ [16J. 

The hopping amplitudes for atomic orbitals can be, in general, complex ^T\. However, 
for optical lattices, coupled waveguides, or superlattices, the hopping amplitude, determined 
by the overlap of adjacent on-site (Gaussian or exponential) ground-state wavefunctions, is 
positive, just as it is for s-wave atomic orbitals [21J. A truly non-Hermitian Hamiltonian 
Hpt may be realized in systems with asymmetrical hopping due to an in-plane field [22J or a 
voltage bias [23]. Its Hermitian counterpart H, with position-dependent, parity symmetric 
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hopping, may be realized in evanescently coupled waveguides where the wavepacket evolution 
and two-particle quantum correlations are exquisitely sensitive to the hopping [20J . 

Note that Eq.(j2|, although dependent upon the underlying Hamiltonian H PT , does not 
uniquely determine the transformation matrix M or the Hermitian matrix H. For simplicity, 
we choose M to be real and mi = 1 which implies, via = mk^t^-kl^k, that = 1. 
The resulting real matrix M commutes individually with the parity- and time-reversal oper- 
ators. Since numerical diagonalization of a Hermitian matrix H is faster and more accurate 
than its non-Hermitian counterpart Hpt, in numerical calculations we use its Hermitian 
counterpart H with entries 

H = -If Iat l 1/2 (S ie i6m +S i -ie N - n \ 

-L-Lmn \' J m,' J N —m\ \ u m,n—l^ ~ u m—l,n z ' J i 

where we recall that O^-n = 6 n . In the following sections, we discuss basic properties of 
such a chain, with focus on the energy spectrum and nature of wavefunctions of H when 
the hopping is not uniform. 

Energy Spectrum and Density of States: We start with numerical results for an N = 500 
site chain with hopping amplitude given by tk = tok a where to sets the hopping-energy 
scale. When a = 0, we have a uniform tight-binding chain, the energy spectrum is given 
by E n = — 2to cos(k n ) where k n = mr/(N + 1) for an open chain, and the density of states 
Po(^) = 0(1 — |x|)/27rfo\/l — x<1 diverges near the band edges x = ±1, where x = E/(2t ) 
and 9(x) is the Heaviside function. The left panel in Fig. [I] shows the cosine-spectrum for 
a = (black solid line), a linear spectrum that is obtained when a = 1 (green dotted line), 
and nonlinear spectra obtained when a = 2 (red dashed line) and a = — 1 (blue dot-dashed 
line). As is expected for a tight-binding model, the energy spectra are symmetric around 
zero [23]. These results show that the energy spectrum of the VT -symmetric chain can be 
widely tuned. We note that when a < 0, the eigenstates near the top and the bottom of the 
energy band are localized at the two ends of the chain. 

The right panel in Fig. [l] shows the corresponding (un-normalized) densities of states 
p a (E). It is clear that the density of states changes dramatically from po(E) (black solid 
line) when a^O. When a = 1, due to the linear spectrum, the density of states is constant. 
It develops a single peak at E — and tapers off to a finite value at the band edges when 
a > 1. In contrast, when a < 0, it develops two symmetrical peaks and vanishes at the 
band edges for iV — > oo. We emphasize that when q ^ 0, the system is not translationally 
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FIG. 1. (color online) (a) Left panel shows energy spectra for the robust PT-symmetric chain, 
Eq.Q, with N = 500 sites and a position-dependent hopping amplitude i& = tok a with a = 
{0,1,2,-1}. The energy is normalized by its maximum value. When a = (black solid line), 
we recover the well-known tight-binding chain dispersion E n /{2to) = —cos(k n ). When a = 1 
(green dotted line), we get a linear spectrum. When a > 1 (red dashed line), the energy spectrum 
develops an inflection point at zero energy. In contrast, when a = — 1 (blue dot-dashed line), the 
energy spectrum is linear at the origin, has a steep slope near the band extrema, and develops two 
symmetrical inflection points, (b) Right panel shows corresponding (un-normalized) densities of 
states p a (E). When a = 2 > 1 (red dashed line), p a =2(E) develops a maximum at zero energy and 
it monotonically decreases to a finite value at the band edges. When a = — 1 < (blue dot-dashed 
line), the density of states shows a two-peak structure. When a = (black solid line), we recover 
the well-known result po{E) that diverges at the band edges. These results show that the energy 
spectrum and density of states are widely tunable through the exponent a. 



invariant and therefore the quantum number n is not associated with the momentum. 

We now focus on a = 1 or equivalently tk = t^k for k = {1, . . . , N— 1}. The band-edges in 
this case are given by ±E Q = ±(N — l)i and the uniform level-spacing is AE = E n+ i — E n = 
2t . It follows from Eq.Q that the recurrence relation satisfied by the coefficients of an 

5 



eigenfunction |^ 7 ) = J2k=i fk\ k ) of H PT is 

- *o [kf k+1 + (N + 1- fc)^_J = K,^. (4) 

It is easy to check that /f = C^ 1 = (JV-l)!/(fc-l)!(jV-A;)! satisfies Eq.Q with eigenvalue 
E G = — (N — l)t . Thus the ground-state wavefunction is \i/)q) = ^2f =1 fk\k) ■ Note that 
for k ~ N/2 ^> 1, Stirling approximation implies that the ground-state wavefunction is 
Gaussian near the center of the chain, fff ~ exp[— (k — N/2) 2 /2N 2 ]. The first excited- 
state wavefunction is given by fl = (N + 1 — 2k) f^. It has energy — (N — 3)t , and Stirling 
approximation shows that in the large iV-limit, it carries over to the wavefunction for the first 
excited state of a simple harmonic oscillator. It is straightfoward, but tedious, to construct 
the higher excited states. We emphasize that for every eigenstate with energy — E < 0, the 
eigenstate with energy +E > is given by its staggered version: f k — > (—l) k fk [24|. 

When a ^ {0, 1} an analytical solution for the eigenvalue spectrum Hp?, or equivalently 
H, is unknown. However, the results in Fig. [T] for a > can be qualitatively understood 
with the simplest example of a non-trivial, symmetric, tridiagonal matrix H with real entries 
{a, b, b, a} above the diagonal. The matrix H is similar to a PT-symmetric Hamiltonian Hpx 
of a 5-site chain with hopping parameters {^1,^2,^3,^4} with a = —\Jt\t4 and b = —^2^3- 
The eigenvalues of such a matrix are given by E n = {±\/a 2 + 2b 2 , ±a, 0}. For a position 
dependent hopping t k = t k a , when a = the slope of the energy spectrum at the band-edge 
is smaller than that at the origin, when a = 1 we get the linear spectrum, and for a > 1, 
the slope of the spectrum at the band-edge is larger than that at the origin. 
Localized and Extended Wavef unctions: The PT-symmetric chain with a position-dependent 
hopping is not translationally invariant, and when a 7^ {0, 1} its eigenfunctions are 
not analytically known. To study the evolution of the spatial extent of a wavefunction 
|^) = Y2iLi fi\i) with increasing system size N, we calculate the inverse participation ratio 
(IPR) [25] 

IPR^(iV) = llM (5) 

(EL l/,l 2 ) 

If the IPR remains finite as iV — > 00, the wavefunction is localized whereas for an extended 
state, IPR(iV) oc 1/N V — > where 7] > 0; for a uniform tight-binding chain, a = 0, the 
IPR(iV) = 3/N for all eigenstates. Note that the IPRs for eigenstates with energies ±E 
are the same. The left-panel in Fig. [2] shows the evolution of the maximum and minimum 
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FIG. 2. (color online) (a) The left-panel shows the minimum and maximum values of inverse- 
participation-ratio (IPR) for an iV-site chain with Hamiltonian H, Eq.Q, as a function of position- 
dependent hopping tk = tok a . When a = (black solid circles), we obtain the analytical result 
IPR(iV) oc 1/N. The a = 1 (blue dot-dashed line) and a = 2 (red dashed line) results show that 
the IPRs decrease monotonically with increasing chain size. These results strongly suggest that 
all eigenstates are extended when a > 0. (b) The right-panel shows the IPR results for a < 0. 
The a = —1/2 (blue solid circles) and a = — 1 (red open squares) results show that the minimum 
IPR is essentially independent of a. The maximum IPR saturates to a nonzero value and indicates 
the presence of localized eigenstates with energies ±E. These results show that when a < 0, the 
system has both extended and localized states. 

values of IPR for a chain with N = 10-5000 as a function of a > 0. Note that since the 
chain size and the IPRs span decades, we use the logarithmic scale in Fig. [2j When a = 
(black solid circles) we obtain the analytical result, IPR(iV) = 3/N. When a > 0, we see 
that both the minimum and maximum IPRs decay with increasing chain size, max IPR 
oc N~ Va and min IPR oc N~ la where < r\ a < 7 a < 1, and both exponents r] a and 7 Q are 
monotonically decreasing functions of a. These results strongly suggest that all eigenstates 
of the Hamiltonian H with position-dependent hopping t' k = t [k(N — k)] a ^ 2 are extended 
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when a > 0. The right panel shows corresponding results for a < 0. The minimum IPR 
is essentially independent of a. On the other hand, in a sharp contrast with the previous 
results, we see that the maximum IPR quickly saturates to a nonzero value and indicates a 
localized state. Thus, when a < 0, the system has both extended and localized eigenfunctions. 
We note that these exponentially localized states, at the two ends of the chain, are essentially 
degenerate in energy; so are their staggered counterparts [21]. Thus, there are at least four 
eigenstates that have the same nonzero IPR. The qualitative difference between a > and 
a < cases can be attributed to the hopping: when a > 0, the hopping amplitude increases 
from ~ toN a / 2 at the two edges to ~ t N a at the center of the chain, whereas when a < 0, 
the hopping amplitude decreases from ~ to/N^ a ^ 2 at the edges to ~ t /N^ at its center. 

A better insight into the number of localized states is provided by the dependence of the 
IPR distribution on the size N of the chain. Figure [3] shows the histogram of IPRs for a = 1 
(left column) and a = — 1 (right column). When a = 1, as N is increased tenfold from 
N = 500 (bottom-left panel) to iV = 5000 (top-left panel), the entire IPR distribution shifts 
to smaller values. In contrast, when a = — 1, even as N is increased tenfold, the IPR values 
for a few (localized edge) states, indicated by the red oval, are unchanged while the IPRs 
for the rest shift to lower values. 

Discussion: PT-symmetric lattices with a uniform hopping and imaginary impurities have 
an extremely fragile PT-symmetric phase [T2l IT3] . In this paper, we have presented a 
tight-binding model with non-Hermitian, PT-symmetric hopping, Eq.Q. We have shown 
that, under very general circumstances, this model is always in the PT -symmetric phase, 
and its Hamiltonian is similar to a Hermitian Hamiltonian with position-dependent nearest- 
neighbor hopping, Eq.(|3]) [20J. These results are unaffected by the presence of ubiquitous, 
on-site, Hermitian disorder since it does not induce PT-symmetry breaking. 

Given the robust nature of the PT-symmetric phase in this chain, we have explored its 
energy spectrum, density of states, nature of eigenfunctions, and their dependence on the 
functional form of the hopping amplitude tk = tok a . We find that when a = 1 the energy 
spectrum is linear and gives rise to a constant density of states. We show that the energy 
spectrum is widely tunable by changing a. We find that when a < the system has both 
localized and extended eigenfunctions in the absence of disorder, whereas when a > 0, all 
eigenfunctions are extended. The effect of a Hermitian on-site disorder, then, is identical to 
that in a regular tight-binding model [T71 [2S] . Thus, the physics of the robust PT-symmetric 
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FIG. 3. (color online) The evolution of the inverse-participation-ratio (IPR) distribution with the 
chain size N. Note the logarithmic scale. The top left and bottom left panels show that when 
a = 1, as N increases, the entire distribution of IPRs shifts to lower values. It suggests that all 
eigenstates of the Hamiltonian H, Eq.Q, are extended in the absence of disorder. The top right 
and bottom right panels show that when a = — 1, as N is increased, although most of the IPRs 
shift to lower values, they saturate to a nonzero value for the states shown in the red oval. These 
eigenstates are localized at the two ends of the chain, and each finite value of the IPR is four-fold 
degenerate. Thus, when a < 0, the system has both extended and localized eigenstates in the 
absence of disorder. 



chain with non-Hermitian hopping is extremely rich. 
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